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ABSTRACT 


Wind energy has become ever more popular in recent decades. With its increasing penetration, requirements for grid 
connection have been established. Among these requirements, the low voltage ride through capability is regarded as the 
most challenging requirement found in wind generators which are integrated into the grid. The induction generators are 
behaving differently for symmetrical and asymmetrical faults. In this the brushless doubly fed induction generator 
(BDFIG) is used because of its lower cost and higher reliability when compared with DFIG. This paper explains about the 
dynamic behaviour of BDFIG under symmetrical and asymmetrical low voltage dips supported by MATLAB simulations is 


provided. 
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INTRODUCTION 


The brushless doubly fed induction generator (BDFIG) is used for wind energy generation, especially offshore. Compared 
to DFIG, it has the advantage of being removal of brush gear and slip rings the maintenance cost of the system might get 


reduced and it has the high reliability [1] [2]. 


During recent decades wind energy generation has adept a rapid development. Wind energy commonly 
acknowledged as the clean and environmental friendly inexhaustible energy source that can reduce the usage of fossil 
fuels. With increasing access of wind power wind generators are required to remain connected and reactive current supply 
to the grid during grid voltage dips. Because of the advantages and the demand of wind energy most of the wind turbines 
are attached to the grid are operate at varying speed. At present the market for the variable speed wind turbines is 
increasing and these are capable of generating high power. And among the many requirements the low voltage ride through 
(LVRT)capability is the main requirement in grid code[3]-[5] which means the generators are remain connected to the grid 
to ride through the low voltage faults. There are two types of grid faults which are symmetrical and asymmetrical faults. 
When a symmetrical fault occus, wind turbines need to perform ride through of zero voltage and further grid voltage 
recovery also inject the reactive current to the grid up to the rated current during fault time. Asymmetrical faults are 
introducing a backward sequence to the generator but the current surge is less severe but these are more likely to happen. 


For such applications the Doubly Fed Induction Generator and Brushless Doubly Fed Induction Generators are used. 
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DFIG is the most widely used generator for wind turbines because it achieves the adjustable speed operation for a 
small rated power converter [6] [7].But the existence of brush gear and slip rings the maintenance cost of the system is 
getting higher. And in DFIG the flux from the stator is exposed directly to the grid so if there is any voltage dips then that 
will causes a sudden loss of the magnetization of the machine. DFIG also very sensitive towards asymmetrical faults [9] 


for which DFIG requires a crowbar protection equipment to take the control of transient over currents [11-13]. 


To overcome the above difficulties we are using the brushless doubly fed induction generator (BDFIG) in wind 
turbines. BDFIG has a larger series leakage reactance hence it experiences less transient current compared to the DFIG so 
as a result it is possible for BDFIG to ride through the low voltage fault without the need of an extra protecting equipment. 


So this paper is explained about the analysis of BDFIG in the symmetrical and also asymmetrical low voltage faults. 


BDFIG OPERATION 


The BDFIG stator has two windings with distinct pole pair numbers in order to prevent direct pairing between the 
windings. The rotor employs a special design enabling it to couple both stator windings. The stator windings are named as 
power winding (PW) and control winding (CW). The PW linked directly to the grid whereas the CW linked to the grid 
through the back to back converter which is variable voltage and variable frequency (VV VF) bidirectional converter that is 


handling a fraction of rated power [1]. 


The BDFIG is normally operated under synchronous mode i.e., in doubly fed mode in which shaft angular 
velocity determined by the excitation frequencies of the PW and CW of the stator, which is independent of the torque 


exerted on the machine, and can be expressed as the 


— @1t@2 (1) 


, PitD2 


where w,and Wz are the excitation angular frequencies supplied to the two stator windings. 


When the CW angular frequency w, equals to zero, the angular velocity of the shaft becomes the natural angular 


velocityw,. 


oy = ace (2) 


P1tp2 


The vector model of the BDFIG aligned in an arbitrary rotating reference frame of w is defined as [17] 











Vi, =Ri +S4+jo 6, (3) 
6 Ah ake (4) 
Vo = Roly +22 + j(@ — (Pi + pz) or) b2 (5) 
b>, =Lal, + Lol, (6) 
V, = Rrly + 224+ j(o — pros), (7) 
fj eel, elisa (8) 


Where R,, Ro, Ly, Lz, L4,,L2,,L, are the PW resistance, CW resistance, PW self-inductance, CW self-inductance, 


mutual inductance between PW and rotor, mutual inductance between CW and rotor, rotor self-inductance respectively? 
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The active power and reactive power of the power winding (PW) of the stator are expressed as 


P, = =Re[V, hy] 


~ 2 
3 — 
Q= 5 im[V11,] 


In the steady state all the vectors are rotate at the speedw,, which is set by the grid frequency because the PW 


connected to the grid. 
V, = Vi let (9) 
Where |V,| is the magnitude of the PW voltage. 


L +LorL41 
— _tbirditberlil (10) 


ri 2 LiRr 
Lilr-Lirt+sjp, ae 


From equations (4),(6),(7) &(8) the CW flux linkage can be represented in terms of PW flux linkage and CW 














current as 
2 2, bil2Rr 
_ lazise LyLaLy—-LyLér 2lMar*Sjp or 7 u 
P22 4 eer Qi + LiL?» ae 2 (11) 
Ir Ss jpyor tor “16 s-jpi or 


Where s is the differential operator ae 


The term s — jp; m, is very large so the reciprocal will be neglected. Then the control winding voltage will be 











LyLgLy—L4L3,-L2L4, . Ligliox : 
Vz = 1,R2.+ ao ne. — jy + p2)o, n+ iP a (S- j (Py + P2)Wr)P1 
1¢r a7 1r7 *1"r 
Vo = Vyo + Eys (12) 
Where 
Vyo= (R2+L)(s — j(P1 + Pz) © Lily (13) 
— _lLirber : 
Bor = gi (Sj Pi + Par) oy (14) 
a2 2 
p= ae as (15) 


From above the CW voltage is split into two terms which are EMF induced by the PW flux it is referred as Ey, 
and Vy,is defined as the voltage drop across the CW resistance and equivalent leakage inductance between the PW and 


CW which is caused by the CW currentl,. 
For simplicity the PW resistance is ignored for analysing the behaviour of the BDFIG under the faults. Hence the 


equation (3) can be written as 


d : 
Vy, = 3 = jor91 (16) 


Substituting (9) and (16) in (14) gives 





Lirlar ©1~-(PitP2) or \V,Jeloxt 
1 


E.,= 
61 [2 -LiLy 4 
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— _lirbar jo4t 
= ts, IVile (17) 





Where Sy is the slip of the BDFIG? 
The EMF induced Eg, also written as when we are taking CW as a reference frame as subscript (2) 


(2) _ _Lirlar —i(@1-@1+P2)or)t 
Boi ae Sn|Vile r (18) 


And the voltage drop due to CW current with respect to CW reference frame is 
Va2) = (Ro+Li(s — jr + P2) o LLY” (19) 


Hence by combining equations (15) & (16) the CW voltage in terms of PW flux and the voltage drop caused by 


the CW current by considering the converter point of view is 
(2) _ pl) (2) 
yO = £2 +6 (20) 


The above equations (15), (16) and (17) shows that normal operation of BDFIG the CW voltage depends on the 
PW voltage and its gain is proportional to the CW frequency which is determined by the rotor angular velocity w,. The slip 


range for the fractionally rated converter is from -0.3 to 0.3. 
BEHAVIOR OF BDFIG UNDER SYMMETRICAL VOLTAGE DIP 


Based on the grid code the voltage dip occurs rapidly and grid faults are common in nature. In this symmetrical short 


circuit of the grid is considered during this the voltage collapse to zero. Consider the fault occurs at the time of t=to, 


Then 


Jwyzt; 
a has a 
, = *0 


CW open circuit 


Whenever the fault occurs the PW flux linkage is confined in the machine and the rotor windings continue to cut 


the trapped flux at the prefault speed and the larger EMF E,,, which is observed at the CW though the rotor of the BDFIG. 


In this we have to assume the CW is in open circuited. Hence the current flow through the winding is zero. The flux 


is a state variable and it cannot be discontinuous, hence the PW flux evolves from initial value to zero exponentially, then 


ae ‘fie it, 
1 
a Val gest ad (22) 
re et if t =to 
By using the equations (5),(3) and (14) the time constant of flux decay is derived as 
_ Laly-Lir 
aad (23) 


By the (22) before the appearance of fault (t<t)) the steady state flux vector rotates at the synchronous angular 
velocity w,with constant magnitude. Whenever the fault appears then the flux vector stops at tgand then it decays 


exponentially after that time period (t>tp ). 


Impact Factor (JCC): 5.6125 NAAS Rating 2.96 


Dynamic Behavior of the Brushless Doubly Fed Induction Generator under Symmetrical and Asymmetrical Faults I5 


Then by substitute (21) and (22) in the (17) for t = tp, with respect to PW reference frame the EMF induced in the 





CW is 
E..= _Lirlor_ (-j(p, + 24 
gi = fn iale calale —j@1 + pz); o4) (24) 
L rL 7 6 —(t-to) 
et (-sn)Vil estee™ = (25) 


By considering the CW is open circuited the CW current is zero then CW voltage will be equal to the induced 
EMF. Because of the flux the CW voltage is a vector fixed to the PW stationary frame where the magnitude decays 
exponentially. Hence by the converter point of view the CW voltage vector reverses the rotation and has the angular 


velocity of (p, + p2)w;. 


@) _ p) 
V, =Eyi 





LirL —(t-to) . 
_ st 2r as sn) Vi él ortoe el (P1+P2) ort (26) 





From the above equation the CW maximum voltage is obtained when the fault obtained is given by 





Wel = a (1-sn)MiI. (27) 


The maximum voltage is directly proportional to (1-s,,). Meanwhile the slip of the BDFIG is usually ranging from 
-0.3 to 0.3, to a speed range of +30% hence the maximum CW voltage when the fault happens and the machine is running 
at the rated speed is 130% of the natural speed. By considering the equation the voltage is directly proportional to the slip 


ls, | then the maximum voltage at the fault is 4.3 times the CW voltage before the fault. 
Analysis of BDFIG under Asymmetrical Voltage Dip 


In this case also we are taking the CW is on open circuited which means the current through the CW is zero (J, = 0). For 
the analysis of asymmetrical faults we are considering the sequence component method [18][19]. Whenever a fault appears 


at t = appears at t = ty then the PW voltage will be defined as, 


Vi leeripe-< ty 
Y= jwyt —jw,t ‘ (28) 
[Vi ple? V4 1V,,|e 7°, if t >to 
Where then subscripts p, n represents the positive and the negative sequence which are decomposed as 
2n V. 
V. = ie)“ 
ee [i ¢ aes | (29) 
I inl Je ae Vic 








Via = Vsin w,t 
. 20 
Where Vip = Vsin( (a) zt + =e) 
Vic = Vsin(w,t+ =) 


The flux generated by the forward sequence is rotating at the synchronous speed and the flux generated by the 
backward sequence is rotating at the synchronous speed but it is in the reverse direction. The equation for the flux of 


positive and negative component is obtained as 
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ipl oj 
Pip = ee (30) 
b af = IVinl ev) oy4t (31) 


-jo4 
The positive and negative sequence components of the PW flux linkage in the steady state can be varied with the 
type of fault which is obtained. In this paper we will discuss about the three types of asymmetrical faults which are phase 
to phase short circuit (p-p), one phase to the ground short circuit (p-n) and phase to phase to ground short circuit (p-p-n). 


The PW of the BDFIG is connected in A then the PW voltage and flux linkage are expressed as 


Phase to phase short circuit (p-p): 





























0 
Via V3 . 311. 
Vip = ae sin(w,t + =) (32) 
Vac By sin(w,t + =) 
1 Veit 
[Pm] = + |? ews (33) 
Pin jo4 1Vei wyt 
2 
For phase to ground short circuit (p-n): 
W3ap xs i 
Vie a sin(w@,t+ 7) 
Viv} =| Vsin(w t+) (34) 
Vic v3 7 tT 
ai sin(@,t+ =a) 
2 jot 
ie eal ae (35) 
Pink jer |tye-KotD 
3 
Phase to phase to ground short circuit (p-p-n): 
" 0 
la ener 3x 
Vie By sin(@ t+ 5) 
1 joyt 
o) ) ai 3ve 1 
= (37) 
o in 


jor 1Ve-i yt 
3 


Whenever the transient state happens the PW flux linkage is continuous from one state to the other state and also a 
transient zero sequence is introduced to connect these two states although there is step change in voltage[20] when the fault 


happens at t = to. So the flux linkage for zero sequence @ 1, is added then the PW flux linkage is expressed as 
Pip t Ointoiz = by (38) 
Substitute the equations (30),(31) in (38), then the zero sequence flux linkage written as PW reference frame is 


IViI-lv a o4tyy ew w to —(t-tg) 
bi, = (Wale Wangs Maile a? e = (39) 


J1 
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Since we are considering the CW is open circuited then the CW current is zero, the CW voltage is the EMF 


induced From CW stationary reference frame from the equations (14),(20),(38) and (39) the voltage is 


Ve? = Bip + Eoin t+ Bois (40) 
Where 
o1p = oe inl Vaplem et (41) 
nee Cea (42) 
Biz = Eo seal pleat (43) 


From converter side the open circuited voltage that contains all the three components during an asymmetrical fault at t > to 


(2) 
o1p 


@) . 


Steady state positive sequence of E4;,: 


In this sequence E*;, rotates at an angular velocity of w,— 


(p; + pz) w ,,with slip of |s,|. It is aligned with the prefault EMF but with a smaller magnitude due to forward sequence of 
grid voltage which is under fault !V,,| which is smaller than the rated voltage IV,|. Usually the BDFIG is designed with a 
speed range of +30%, hence the slip is -0.3 to 0.3. 


(2). 


gin: in this sequence EC ,fotates at an angular velocity of —w,— 


Steady state negative sequence of E 
(p; + pz) © ,,with slip of |2 — s,|-In this, the scaling factor is much larger than the positive sequence and maximum is the 
2.3 when the BDFIG is running at high speed. That means the speed is 130 % of natural speed. The negative sequence 


component voltage IV,,,| is smaller than the grid rated voltage. 


Transient zero sequence of E ae In this sequence E ee ,rotates at an angular velocity of —(p, + p2)w,,with slip of 
|1 —s,,|. In this case also the value is more than the positive sequence with a maximum value of 1.3. The magnitude of 


zero sequence of the CW voltage is obtained by the time whenever the fault occur and decays with the time constant of T 




















_72 
which can be obtained as “27-2 
RyLy 
Table 1: Prototype Machine Specifications 
Parameter- Value Parameter- Value 
Machine rating-250KW L, -0.105 H 
PW/CW pole pairs-2/4 Lz -0.382 H 
PW rated voltage-690V(50HZ) L,, -0.004 H 
PW rated current-178A(line) L5, -0.006 H 
CW rated voltage-620V(50HZ) L, -2.602x 10-*H 
CW rated current-84A(line) R, -0.079 ohm 
Operating speed-350~650r/m R,-0.621 ohm 
Rated torque-3670Nm R, -1.770x 107*ohm 
EXPERIMENTAL RESULTS 


The results of the simulation which are constructed under the MATLAB/SIMULINK to verify the mathematical model of 
the brushless doubly fed induction generator under the symmetrical and asymmetrical low voltage faults. A 250KW rating 
BDFIG is taken for this analysis. When the fault happens at the time t=0 sec then the EMF induced in the control winding 


of the BDFIG is raised and after fault clearing it is decaying exponentially. The Figure shows the behaviour of the BDFIG 
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under symmetrical fault the control winding EMF and flux waveforms. The Figure 2 shows the behaviour of the BDFIG 


under the asymmetrical fault which is one phase to ground fault the EMF induced in the control winding, EMF induced in 


the rotor winding of the machine. In figure 3 shows the behaviour of BDFIG under asymmetrical phase to phase fault in 


this also the EMF of the rotor winding is decreased exponentially after the fault clearing. 


Impact Factor (JCC): 5.6125 


Source Voage (volts) 
eo 8 


um aA mE Z iw 
foo 
— a 




























Hi 





i vi 
i 





mm | 


i 
me 


hi 
i 


m 
Kl 


Nh 


i 
t i 


i 
“ 


i i 
i 


i i 


ae 
ti wl 


i 


phi : 


Hh | 


i 

















Hh ! 


























Halil 
















































Boch Wit i i 


sien 
wi 


nh 

















Figure 1: Symmetrical Fault. 
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Figure 2: Asymmetrical Fault (Phase To Ground). 
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Figure 3 Asymmetrical Fault Phase To Phase. 
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CONCLUSIONS 


This paper analysed the behaviour of the BDFIG under symmetrical and asymmetrical low voltage faults. The 


asymmetrical faults that consider here are phase to ground and phase to phase low voltage faults. The experimental results 


proved that the without any need of external hardware such as crowbar circuits BDFIG is able to meet the grid regulations. 
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